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Three-dimensional  optimal  missile  avoidance  Is 
analysed  with  a linearized  kinematic  model.  The 
solution  requires  maximum  load  factor  and  the  prob- 
lem is  reduced  to  optimal  roll  position  control 
having  two  phases:  (1)  orientation  oi  the  lift  vec- 
tor into  the  optimal  evasion  plane,  (2)  rapid  180° 
toll  meanuvers  governed  by  a switch  function.  For 
circular  missile  vectograms  the  plane  of  optimal 
evasion  is  perpendicular  to  the  line  of  sight.  Eva- 
ding from  roil  stabilized  missiles  of  rectangular 
vectogram,  further  advantage  can  be  taken  maximizing 
the  target -missile  maneuver  ratio.  Bounded  roll- 
rate  reduces  the  miss  distance  but  does  not  affect 
the  optimal  evasive  maneuver  structure. 


Nomenclature  > 

lateral  acceleration. 

system  matrix  (14). 

single  channel  matrix  (15). 

control  matrix  (18). 

variational  Hamiltonian. 

pay-off  function  (48). 

constant  of  true  proportional 

navigation  (9). 

miss  distance  (46). 

effective  prop.  nav.  ratio  (9). 

roll  rate  (control  variable). 

relative  distance 

switch  functions  (55),  (68). 

time. 

nominal  time  of  flight  (8). 
control  vector, 
ve  loc  i t y . 

state  vector  components, 
components  of  K,  perpendicular  to  the 
initial  1 i ne  of  s i gh t (12). 
dynamic  similarity  parameter  (92). 
costate  dependent  coefficients  (52), (53). 
normalized  time -to- go  (63). 
costate  vector  components, 
missile-target  maneuver  ratio  (30). 

H M ••  •«  n 

or  a 

single  channel  (33). 
missile  time  constant, 
roll  angle. 

normalized  target  roll  rate  limit, 
azimut  angle, 
angular  velocity. 


Super sc r i pts 


3-1)  vector. 

nondimensional  variables, 
transposed  of  a matrix, 
optimal  control  functions, 
time  derivative. 


Subscripts 
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commanded  value, 
final  value, 
i ndex. 
missile . 
required  value, 
line  of  sight . 
target . 

initial  value, 
column  vector. 


L _ Introduction 

The  missile-aircraft  pursuit -evasion  problem  can 
be  either  formulated  as  a zero-sum  differential  game, 
or  decomposed  to  two  reciprocal  optimal  control 
problems  whose  respective  objectives  are  to  deter- 
mine: (a)  optimal  guidance  laws  against  maneuvering 
targets,  (b)  optimal  evasive  maneuvers  from  guided 
missi 1 es . 

Regardless  of  the  formulation,  the  problem  is  of 
an  inherent  complexity.  The  relative  pursuer-evader 
kinematics  are  expressed  by  a nonlinear  three- 
dimensional  vector  equation.  Botli  vehicles'  dynamics 
are  expressed  by  sets  of  nonlinear  differential  equa- 
tions. Moreover  the  guidance  law  of  the  pursuer  is 
implemented  by  a rather  complicated  transfer  func- 
tion. Hie  exact  solution  for  each  of  the  alterna- 
tive formulations  requires  the  solution  of  a non- 
linear two-point  boundary  value  problem  of  very  high 
dimension.  Hie  computation  of  such  a solution,  al- 
though feasible.  Is  so  time-consuming  that  it  makes 
this  approach  impractical  for  systematic  studies. 

For  a systematic  analysis  which  is  necessary  to 
create  an  insight  into  this  complex  problem,  simpli- 
fied analytic  solutions  are  required.  Important 
simplification  can  be  achieved  by:  (a)  neglecting 
guidance  dynamics,  (b)  restricting  the  motion  in  a 
plane,  (c)  trajectory  linearization. 

It  turns  out  that  the  attractive  assumption, 
made  by  neglecting  the  dynamics  of  the  pursuer, 
yields  seriously  misleading  resul ts 1 * 2 • * . As  a con- 
sequence of  this  assumption,  the  direction  of  the 
optimal  evasive  maneuver  is  constant  and  is  deter- 
mined by  the  initial  or  terminal  conditions.  More- 
over, if  the  pursuer's  maneuverability  is  sufficient, 
the  final  miss  distance  is  always  zero*. 

Most  analytic  studies  in  the  past  used  two- 
dimensional  models'4"*0.  Whenever  guidance  dynamics 
were  considered  (even  if  by  an  approximation  of  a 
first  order  time  constant  or  a pure  time  delay)  an 
oscillating  or  "bang-bang"  structure  of  the  optimal 
evasive  maneuver  became  apparent.  It  was  also  shown 
that  optimal  evasion  can  guarantee  non-zero  miss 
distance  even  from  a pursuer  of  unlimited  maneuvera- 
bility7 or  from  one  of  an  optimal  guidance 
strategy^.  It  has  been  indicated  however10  that  the 
validity  of  2-D  analysis  is  limited  to  near  "head- 
on"  or  "tail-chase"  engagements.  For  other  initial 
conditions  three-dimensional  analysis  is  required. 
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"bang-bang"  structure  of  the  optimal  evasive  maneu- 
vers, trajectory  linearization  is  a good  approxima- 
tion for  a wide  range  of  parameters. 

Hie  objective  of  this  paper,  imtivated  by  the 
above  mentioned  results,  is  to  analyse  the  problem 
of  optimal  missile  avoidance  using  a three- 
dimensional  linearized  kinematic  model.  Analysis  is 
based  on  the  following  set  of  assumptions: 

(1)  Pursuer  and  evader  are  both  considered  as  con- 
stant speed  mass  points.  (2)  The  pursuer  is  a 
homing  missile  launched  against  an  initially  non- 
maneuvering evader  (target)  in  a collision  course. 

(3)  Relative  pursuer-evader  trajectory  can  be 
linearized  around  the  initial  line  of  sight. 

(4)  Pursuer  and  evader  both  have  perfect  information 
on  the  relative  state.  (5)  Gravity  can  be  neglected 
for  both  vehicles  (not  effecting  relative  trajec- 
tory). (b)  Hie  pursuing  missile  has  two  identical 
and  Independent  guidance  channels  to  execute  propor- 
tional navigation  in  two  perpendicular  directions  in 
a plane  normal  to  the  line  of  sight  (true  propor- 
tional navigation11).  (7)  Hie  dynamics  of  each 
guidance  channel  is  assumed  to  be  (for  sake  of  sim- 
plicity only)  of  first  order.  Hie  validity  of  the 
first  five  assumptions  and  the  effects  of  more  com- 
plex pursuer  dynamics  are  discussed  in  detail  in 
Ref.  10. 

Based  on  l he  above  listed  hypotheses  the  3-D  mis- 
sile avoidance  is  formulated  as  a fixed  duration  op- 
timal control  problem  maximizing  a terminal  pay-off 
(the  square  of  the  miss  distance).  The  control  vari- 
able is  the  lateral  acceleration  vector  of  the  eva- 
ding airplane.  Hiis  acceleration  is  perpendicular 
to  the  velocity  vector,  its  magnitude  is  bounded  by 
the  limit  load  factor  (or  maximum  lift)  and  its 
direction  is  controlled  by  the  airplane's  roll- 
orientation. 

First  a mathematical  model  of  unbounded  missile 
maneuverability  and  infinite  airplane  roll-rate  is 
used.  This  linear  formulation  leads  to  a closed 
form  solution  and  provides  the  basic  insight  into 
the  problem.  In  consecutive  steps  saturation  of 
missile  acceleration  and  realistic  roll  dynamics  of 
the  evading  airplane  are  introduced. 

Hie  solutions  obtained  by  the  linearized  3-D 
model  are  compared  both  to  the  prediction  of  a 2-D 
linearized  analysis10  and  to  results  of  complete 
non-linear  (6  degrees  of  freedom)  simulation. 

II  Mathematical  Modelling 
Three-dimensional  vector  formulation 


A three-dimensional  pursuit-evasion  is  described 
by  the  vector  equations 


* - *r  - 

a x i 


i5k 


I5I3 


(1) 

(2) 


The  acceleration  command  of  the  pursuing  missile  is 
given  by  Assumption  6 as 

(&,  x A) 

•c  N Itl 


frt) 


(3) 


while  the  actual  acceleration  is  determined  (see 
Ass.  7)  by 


tVm  ♦ VM 


ftd. 


(4) 


The  acceleration  of  the  constant  speed  evader 
(the  target)  is  normal  to  its  vector  velocity 


Vj  ■ (&j«  x ^) 


(5) 


Trajectory  linearization  around  the  initial  colli- 
sion course  (Ass.  263)  yields 

\l\  ^ - VR  a const  (6) 

and  as  a consequence 

|R(t) l-lRol-  VRt  - v„(tf-t)  (7) 

determining  the  final  time  of  the  pursuit  by 

tf  - Ifjl  (8) 

Substituting  (6)  and  (7)  Into  (3)  and  defining 

KN  A.  N'Vr  (9) 

yields 

($m)c-  (7"ty?  (A(t)+(tf-t)A(t))  do) 

The  system  of  differential  equations  (1),(4),  (S) 
and  the  linearized  feedback  relation  (^0)  determine 
£he  9 components  of  the  vectors  R(t),  v T ( t ) and 
VM(t),  if  initial  conditions  and  the  target  angular 
velocity  vector  0j(t)  are  given.  In  the  problem  of 
optimal  missile  avoidance  this  last  quantity  Is  the 
control  variable. 


Non-dimensional  scalar  equations  - linear  case 

The  Initial  collision  plane  (Ass.  2)  Is  taken  as 
glaije  of  reference  for  the  direction  of  the  vectors 
R,  Vh  and  y (see  Fig.  1). 


Fig.  1.  Initial  collision  geometry 

By  choosing  the  X axis  of  the  coordinate  aystem 
to  coincide  with  the  Initial  line  of  sight,  only 
those  components  of  the  relative  motion  which  are 
normal  to  this  direction  have  to  be  considered.  The 
linearized  equation  of  motion  along  this  axis  Is  al- 
ready solved  by  (7).  The  state  vector  la  reduced  to 
be  of  six  components 

XT-  (xi...  Xe)^  (y,y,yM,z,z,zM)  (11) 

"y"  and  "z"  being  the  relative  displacements 
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(17) 


perpendicular  to  the  initial  line  of  sight  (aee 

71*.  2) 


VT  “ yM 


(12) 


Fig.  2.  3-D  pursuit-evasion  geometry 
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Fig.  3.  Roll  orientation  of  missile 
maneuver  planes. 

Moreover  if  the  Y and  Z axes  of  the  coordinate  sys- 
tem are  oriented  to  be  in  the  missile  maneuver 
planes  (see  Fig.  3)  the  state  equation 


x - A( t) x + B u 

has  a decoupled  matrix  A(t) 
f A | (t)  ’ 0 


with 


A(t) 


Ai(t). 


0 

0 

N' 


V - - 

: A,(t>j 

i 
0 
N’ 


T(tf-t f x(tf-t) 

The  control  vector  u has  two  components 

uT  - (ui,u2)  - (ajainfp  a^coa^)  (16) 

a^  being  the  evading  target  lateral  acceleration 
constrained  by 


0 f |aT|±(aT)B 


and  the  roll  angle  is  Measured  relative  to  the 
collision  plane  and  B la  given  by 


0 , - coaxToCoe»M  ,0,0,  cosx^sin^  , 0 
0 , ain$M  ,0,0,  coa^  , 0 


(18) 


Use  of  nondlmenslonal  variables  reducea  the  number 
of  indeoender 
results'2.  1 
distance  by 


of  independent  parameters  and  yields  generalized 
results'2.  Introducing  nondlmenslonal  time  and 


t - t/T 

ft-  R 


(19) 

(20) 


leads  to  normalize  velocity  components  by  1(07)^,, 
and  accelerations  by  (sf)max  • As  a result  the 
state  equation  (13)  Is  transformed  to 


dft/dt  - ft(t)x  + B u 

with  a nondlmenslonallzed  state  vector 


vr 


•v  d£  d2yM  ^ dz 

y*  d?  • ■ z’  a • 


(21) 


(22) 


and  a normalized  control  vector 

3(1-  (aTs  1 n$T , syiospj)  (23) 

The  decoupled  structure  of  the  state  matrix  is 
preserved  with 


Ai(t)  - 


0 

- 1 
-1 


(24) 


(13) 


(14) 


(15) 


[<*rb  (trb 

Nonlinear  effects 


The  state  equation  (13)  or  (21),  describing  sys- 
tem dynamics,  is  linear  due  to  the  Implicit  assump- 
tions of  unlimited  missile  maneuverability  and  In- 
finite target  roll-rate.  A more  realistic  mathema- 
tical model  has  to  consider  the  constraints  on 
these  variables.  The  state  equation  Including  such 
effects  will  no  longer  be  linear. 

Limited  missile  maneuverability 

When  missile  maneuverability  constraints  are 
taken  Into  account  It  Is  necessary  to  redefine  the 
state  vector  and  the  state  equation.  The  compo- 
nents of  the  lateral  acceleration,  y*  and  z^,  are 
to  be  replaced  by  their  required  value  (y|()E  and 
(zM)r  In  nondlmenslonal  form 


*3 


(®T^ max 

<«H>r 

(*T^max 


(25) 

(26) 


As  these  variables  are  trot  affected  by  the  con- 
straint, the  differential  equatlona  for 
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dxj/dt  and  uKfa/dt  remain  unchanged . 

The  relation  between  the  components  of  the  accel- 
eration, which  are  subject  to  constraints  and  their 
required  value  can  be  expressed  by  the  nonlinear 
saturation  function  defined  as 


|slgn|~ ) If  |a|>|b| 

Consequently  the  stale  equation  will  not  have  the 
linear  form  of  (21)  but  has  tu  be  written  as 

dx/dt  - £ (x  , t)  + B u (28) 

For  such  saturation  two  alternative  formulations 
exist  expressed  by  two  different  vectograms: 

a.  Circular  (Isotropic)  vectogram  (see  Fig.  A) 
showing  that  the  constraint  of  maneuverability  ap- 
plies to  the  resultant  lateral  acceleration 

, (yM)*  +(*„)’ 

•»“  'm*  V'  (aM>max  aat'~U^ 1 (29) 

" max 


with  similar  relation  for  Z|{.  Vectogram  of  this 
type  represents  a roll  stabilized  cruciform  missile 
with  known  roll  orientation. 


YT" ( QM)r 


Fig.  4.  Circular  vectogram  of  missile 
acceleration. 

Such  vectogram  represents  the  maneuverability  of  a 
thrust  vector  controlled  (T.V.C.)  missile  or  of  a 
cruciform  configuration  with  unknown  roll  orienta- 
tion. In  this  model  it  is  assumed  that  saturation 
of  both  guidance  channels  takes  place  simultaneous- 
ly. 

By  introducing  the  miss lie- tar get  maneuver  ratio 
which  is  one  of  the  similarity  parameters  of  the 
problem12,  as 

A ^aM^max 

M - 7— -x (30) 

'aT'max 

(29)  can  be  written  using  (21))  and  (26)  in  a non- 
dimensional  form 

r (square) 


b.  Rectangular  (square)  vectogram  (see  Fig.  5) 
Indicating  that  saturation  may  occur  In  each  gui- 
dance channel  separately: 


(yM>r 

' <5^  max  8at  (yj 

•'Umax 


Z 

(ZM)max 

1 

1 

* X 

Fig.  3.  Vectogram  of  cruciform  missile. 

Defining  the  relative  maneuverability  of  the 
guidance  channels  as 

„ A (yM>max  (*M>mex  .... 

’ ^®T^max  (“T^msx 

enables  to  write  (32)  in  nondlmenslonal  form 

d2y  % 

^7  “ Hi  sat  {^M  (34) 

dt 

A similar  expression  holds  for  the  second  channel. 

For  rectangular  vectogram  the  actual  missile 
target  maneuver  ratio  p depends  on  the  orientation 
of  the  target  acceleration.  It  can  be  easily  seen 
that 

Pi  I v I ^ hi  (35) 

Limited  target  roll-rate 

Realistic  description  of  the  evader's  roll  dyna- 
mics requires  thst  its  roll  orientation  +-  be  con- 
sidered not  as  a control  but  an  additional  state 
variable . 

*7  “ 9T  (36) 

The  roll  dynamics  can  be  expressed  by  several 
alternative  formulations: 

1)  The  control  variable  Is  the  target's  roll- 

rate 

i7  - PT  (37) 

subject  to  the  constraint 

|PTI  II  (38) 

2)  Closed  loop  roll  control  based  on  the 
required  roll  orientation 

(x7)r  - k^  l(*7)r  - X7J  (39) 

subject  to  saturation 


*7  ” ^T^max  8at  ^(I-)  I 
Tmx 


(32) 
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B/*,  »,  Vj  it  • [£({)£  * t £)  CX/t 

' an  0*  mfiUD,  **fi*r*»  lo|r  th*  part  loBapaoBcot  of 
tii*  control  yarlaol**,  a* 

B * ...  t>  ♦ ij  [t  ,*1i4t+.  jeo**,;  (51; 

with 

<j«  (//*f^Oi^lM|f->2  coM^/  (52> 

tj.»  I>7»1i4k*'i  coa«||J  (53) 

Tli*  optlaai  rootrol  variable*  a*  and  ♦*  have  to 
*Mi/Jml£*  tli*  Maid  1 Ionian,  yielding 

(a'7l*  - | lalgn  S,  ♦ 1)  (54) 

wl  tli 

Sj  ^ tj,i0g7  * <^oo*d7  (55) 

«r*d 

( ♦■) ) * - tg-Mii/^J  (56) 

b'jiiMi  ltutlng  (56j  Into  (55)  lead*  to 

Sj  - k ♦ «22)  > 0 (57) 

k belriK  « ponftfve  constant  of  proportionality 
deternlrtf'd  by 

k*  - (4,*  -^j2)-1  (58) 


fafuMt  i o fin  and  07)  Indicate  directly  that  for 

optlnaJ  aibbile  avoidance  naxltuoa  lateral  load  fac* 
tor  tiaa  to  be  alwaya  uaed. 

'/lie  component*  of  the  coatate  vector  ^ Involved 
In  Sj,  4 j and  are  determined  by  the  adjoint 


e/juat  Ion 

Aj  A - - bH/^i  (59) 

wlfh  the  te mi  rial  conditions 

A<V  - - <«»> 

reaulllnK  in 

- - 2'f(tf)4  - 2 yf 

*%(t|)  - - 2»(tf)A-  2»f  (*u 

»,('(■,)  - o I - 2, 3, 5, 6 

(59)  ylrlilH  for  a linear  «y«ten  aa  ( 21  > 

dl/d t - - ^('<T)  A (82) 

Int  roduclnK  the  rairmallaed  tlme-to-go 

« ^ - t (83) 

leada  to  traiiMform  (60)  and  (62)  to 

d^/dn-  Wif))  ^ (86) 

with  the  Initial  condition* 

Ao  - - M|6-o  (65> 
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The  system  of  equations  (64)  can  be  reduced  to  two 
Identical  scalar  differential  equations  of  the  form 

O ♦ +«'  5^  ■ 0 i“3,6  <“> 

wiiich  were  solved  in  a closed  form  in  Ref.  10. 

'tills  solution  combined  with  (65)  yields 

Aj(e>  - Aloe'6  [1  - Pj(0)J 

Xl+l<®>“Alo0e_8  I1  - pl+l<e>l  (67) 

*l+^0)-*lo82e'6  II  - Pi+2<8>1  w>4 

P((6)  are  functions  of  0 depending  on  N'  only. 
For  Integer  values  of  N'  these  functions  are  poli- 
noaials  of  the  order  (N'-2)  (see  Table  A-l  In 
Ref.  10). 

As  a consequence  of  (67)  and  (61)  the  time 
dependent  part  of  A;  and  As  are  Identical! 

*\i  ly,  ^ 

*2(t)  * - 2yfS2(t ft) 

n.  a.  „ * (88) 

As(t)  - - 2zfS2(tf-t) 

Substituting  (68)  into  (52)  and  (53)  yields 

r\j  r\j  i\j  <\j 

2cosX'i0lyfCos(|>M-  z f s in4> je| ] S2 ( t f ”t)  (69) 

6?*  - 2[yf8in4>j,(4-  z^cos^^lS? (tf-t)  (70) 

The  expressions  in  the  brackets  are  the  compo- 
nents of  the  miss  distance  vector:  the  first  in  the 
plane  of  the  initial  collision  and  the  second  in 
the  direction  perpendicular  to  this  plane. 

As  by  (54)  and  (57) 


maximizing  the  Hamiltonian  is  equivalent  to  maxi- 
mizing Sj.  Inspection  of  eq  (57), (69)  and  (70) 
makes  it  obvious  that  for  any  given  miss  distance, 
maximization  is  achieved  only  by 


and  as  a consequence 


This  clearly  indicates  that  the  direction  of  the 
optimal  evasive  maneuver  has  to  be  perpendicular  to 
plane  of  initial  collision.  The  optimal  roll  orien- 
tation, either  <£7^  or  $7*11 

cos(<frT)*(t)  - - sign[S2(t f-t) ] (74) 

is  determined  by  the  sign  of  the  switch  function 
S?  which  is  identical  to  the  one  obtained  in  2-D 
analysis10. 

Case  of  circular  missile  vectogram 


nonlinearity  and  the  optimal  control  functions  are 
given  in  this  case  also  by  (52)  and  (53) • Moreover, 
as  saturation  takes  place  in  both  guidance  channels 
simultaneously,  the  time  dependent  parts  of  A2  an<J 
A5  remain  identical  (although  they  are  different 
from  the  costate  variables  of  the  linear  case) • As 
a consequence  both  (71)  and  (74)  hold,  yielding  the 
same  type  of  "bang-bang"  maneuver  perpendicular  to 
the  initial  collision  plane  as  for  unlimited  missile 
maneuverability . Hie  switch  function  governing  this 
maneuver  is  however  different  from  the  one  obtained 
in  closed  form  for  the  linear  problem. 

When  missile  saturation  takes  place  in  the  ter- 
minal phase  of  the  pursuit  (it  was  shown13  that  it 
always  occurs  in  this  phase),  the  state  equation  of 
the  system,  and  as  a consequence  the  adjoint  equa- 
tion, both  change.  As  both  components  of  missile 
acceleration  are  constants  in  the  saturated  phase, 
the  submatrix  in  (24)  is  modi  fed.  For 

Its  second  line  will  contain  only  zeroes.  The  ad- 
joint equations  are  given  in  this  case  as  a function 
of  normalized  time-to-go  0,  as  follows: 

dAj/de  - N’/e2  Ai+2  (75) 

dA^/de-Aj  + N’/e  A1+2  (76) 

dAi+2/de-  - A i+2  (77) 

1-1.4 

with  the  Initial  conditions  (61) . 

From  (77)  and  (61)  it  la  obvious  that  for  6^6, 

A3  and  Ag  are  both  zero.  This  fact  confirms  that  in 
the  saturated  phase  the  required  accelerations  have 
no  influence  on  the  solution.  Consequently 

A1(6<0S)  - (Ai>o  “ const  (78) 

Ai+j(e<0a)  » (A±)o  - 0 (79) 

It  Is  easy  to  see  that,  due  to  the  monotonity  of 
A 2 and  A5,  no  switch  can  occur  when  the  missile  la 
saturated.  It  Is  also  Important  to  note  that,  as 
F(x,t)  in  (28)  has  no  discontinuity  when  saturation 
occurs,  the  costate  variables  remain  continuous  at 
0 - 09.14 

The  time  of  saturation  Is  one  of  the  unknowns 
and  has  to  be  determined  with  the  complete  solution 
of  the  costate  variables,  by  colvlng  a two-point 
boundary  value  problem.  Fortunately,  due  to  the 
"bang-bang"  type  solution  a simple  and  efficient 
search  technique  developed  In  Ref.  10  can  be  used  aa 
an  alternative . The  results  obtained  are  Identical 
with  those  of  a previous  2-D  analysis19.  They  show 
that  the  optimal  switch  function  and  the  resulting 
miss  distance  both  depend  on  the  missile-target  ma- 
neuver ratio  m.  The  dependence  of  the  normalized 
miss  distance  can  be  expressed  approximately  as 

m*(N\  tf,u)  a m*(N* , tf,  -)  + b(N' , tf)  !v2 


Case  of  rectangular  missile  vectogram 


If  missile  maneuver  constraints  are  taken  In 
consideration  the  linear  state  equation  (21)  is  re- 
placed with  one  of  a nonlinear  form  of  (28).  Never- 
theless, this  nonlinearity  does  not  alter  the  struc- 
ture of  the  optimal  solution.  As  the  controls 
appear  separately  in  (28),  the  Hamiltonian  preserves 
Its  separated  form  of  (51).  The  control  dependent 
part  la  not  affected  by  the  saturation  type 


In  this  case,  representing  a roll  stabilized 
cruciform  missile  with  known  roll  orientation,  sa- 
turation occurs  In  each  guidance  channel  separately. 
The  optimal  control  functions  are  determined,  as  for 
the  previous  case,  by  (52)  and  (53).  As  (57)  la  not 
influenced  by  the  saturation,  (71)  remains  valid. 
However,  the  time  dependent  parts  of  A2  and  A5  are 
no  more  identical  and  the  optimal  target  roll 
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orientation  given  by  (56) 

r.A  >*.  Ai  . ™»»'rJA5alnyA2co»*M] 

® T 62  Apsin^H+Ascos^M 

cannot  be  determined  by  simple  inspection. 

It  is,  however,  intuitively  obvious,  in  view  of 
(80),  that  for  maximum  miss  distance  the  effective 
missile-target  maneuver  ratio  has  to  be  ninlmum. 

This  problem  has  a straightforward  geometric  solu- 
tion shown  in  Fig.  6.  it  yields,  for  0<^<ir/4, 
which  is  the  relevant  one  for  cruciform  missiles, 

tg(«£)'  - cosXto  tg$M  (82) 

This  is  a suboptimal  solution,  which  maximizes  the 
projection  of  target  acceleration  on  the  more  sus- 
ceptible guidance  channel,  but  it  can  be  easily 
imple  Dented.  Moreover  It  takes  a definite  advantage 
of  the  known  missile  roll  orientation  by  providing 
always 

1‘eff  i ^2~  Pj  (83) 


t y h ) ma  11 


, (a’i)"  (Zm)hio» 


UQt  )rno,  c°sXr 


/ / effective  target 

TlaOmax  vectogram 


Minimization  of  the  effective  missile- 
target  maneuver  ratio. 


Case  of  li mlted  target  roll-rate 

The  "bang-bang”  solution  obtained  in  the  previous- 
ly discussed  cases  assumes  an  infinite  roll-rate  of 
the  evading  aircraft.  Whenever  the  real  limitation 
on  the  target  roll-rate  is  taken  in  account  the  roll 
orientation  becomes  an  additional  state  variable 
of  the  problem  as  indicated  by  (36).  As  a conse- 
quence the  system  equation  becomes  nonlinear  even  in 
the  absence  of  missile  saturation.  For  such  case 
the  optimal  control  problem  is  formulated  as  follows: 

Given  the  dynamic  system  with  the  state  vector 

•V*  '\v  , 

* = col  |xj . . .x7]=coJ  [x  , (>T1  (84) 

The  state  equation  has  the  form  of  (45)  . its  first 
six  components  are  identical  to  ( 22)  and  the  last 
equation  is  ^ ^ % 

dxy/dt  = PT  (85) 

The  system  is  cont rolled  by 

uT  - 1 Pf]  (86) 

'W 

The  initial  conditions  xQ  are  given  and  the  terminal 
state  is  not  specified . 

'V* 

Find  the  optimal  control  u (t)  subject  to  the 
constraints  (47)  and  (38)  that  maximizes  the  termi- 
nal payoff  (48)  for  the  fixed  hf  given  by  (49). 


The  variational  Hamiltonian  of  the  problem  is 
H’-ATf  - »0<*.  t>-*£x  Si(^,*2,*5>-W,tV7  (87) 

with  S]  written  explicitly  in  (55). 

The  first  optimal  control  variable  a£  is  given 
as  previously,  by  (71).  The  second  control  compo- 
nent maximizing  the  Hamiltonian  haa  to  be  (if  A 7^)) 

(PT)*  - sign  x7  (88) 

The  time  derivative  of  the  new  costate  variable  Is 

dX7/dt  - - 3h’/3x7  - - aT  3Sj/3x7  (89) 

yielding  (as  a function  of  the  noraellxed  tine-to-go) 

dA7/d0  - aT[Si(X2.*5)cosx7-6?(A2,Ai),inx./J  (90) 

with  the  Initial  condition  (A7)0«0. 

A singular  control  is  possible  if  A7-dA7/d6^), 
requiring  by  (90) 

tg  *7  “ t8*T  “ 4(Af.A^)  (’1) 

t Assuming  that  (91)  holds,  the  singular  value  of 
P„  can  be  obtained  from  the  second  derivative. 

This  value  turns  out  to  be  zero.  Comparing  (91)  and 
(56)  indicates  that  the  required  roll  orientation, 
predicted  by  (36)  under  the  assumption  of  an  infi- 
nite roll-rate, does  not  change  by  the  introduction 
of  the  roll  rate  constraint. 

Such  steady  state  (Px«0,  A7m,  tg  Si/4j)  is 
however  very  unlikely.  The  mutual  relations  (85), 
(88)  and  (89),  shown  In  the  block  diagram  of  Fig.  7, 
predict  limit-cycle  type  oscillations  around  the 
equilibrium  value  of  (91).  Thc-e  oscillations  may 
damp  out  if  higher  order  roll  dynamics  are  intro- 
duced in  the  model.  Roll  oscillation  of  small  ampli- 
tude have  no  appreciable  effect  on  the  solution. 

The  major  effect  is  the  reduction  in  the  optimal 
miss  distance  as  the  value  of  the  normalized  maximum 
roll-rate  decreases.  This  phenomenon  was  already 
predicted  by  the  2-D  analysis1®. 


p i*t  m 

(♦r)ma*  — - — r~ 


— K-A  H[x2(t),X5(t)] 


Fig.  7.  Roll  rate  aa  control  variable. 

V Concluding  Remarks 

In  this  paper  the  problem  of  the  3-D  optimal  mis- 
sile avoidance  is  analysed  in  nondime nalonal  form 
for  realistic  missile  and  aircraft  models  using 
linearized  kinematics.  The  solution,  derived  by 
rigorous  mathematical  treatment,  la  presented  in 
simple  geometric  terse , providing  a clear  insight 
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parameter"  introduced  for  nonlinear  kinematics  in  a 
recent  report12  and  defined  as 


« 


% 


Into  this  inherently  complex  problem. 

First,  It  is  shown  that  the  optimul  evasion  does 
not  take  place  in  the  initial  collision  plane.  Ilius 
the  effort  in  J-U  analysis  is  justified.  Neverthe- 
less, the  optimal  evasive  maneuver  is  confined  to  a 
plane  which,  for  circular  missile  vectrogram,  is 
perpendicular  to  the  Initial  plane  of  collision. 
Evading  from  a roll  stabilized  cruciform  missile, 
represented  by  a rectangular  vectogram,  further  ad- 
vantage can  be  taken  by  choosing  a maneuver  plane 
which  minimizes  the  missile-target  maneuver  ratio. 

The  solution  oi  the  optimal  control  problem, 
maximizing  the  miss  distance.  Is  a "bang-bang"  type 
maneuver  witli  the  continuous  use  ot  maximum  load 
factor  of  the  evading  airplane.  It  can  be  therefore 
reduced  to  an  optimal  rol  I -|h>s  1 1 ion  control  problem 
of  two  consecutive  phases:  (1)  Orienting  the  air- 
plane lateral  acceleration  vector  into  the  plane  of 
optimal  evasion.  (2)  ('hanging  the  direction  of  this 
acceleration,  which  has  to  be  maximal,  by  rapid  roll 
maneuvers  of  180°  in  accordance  with  an  optimal 
switch  i unction. 


M * 

3-D  Linearized 

1 5 

2-D  non -linear 

N'  = 4 0 

1 0 

Y,  = 0 13 

If  = 12  5 

4 =0.4 

0 5 

0T  = 4 2 

i. i 

0 15  "30  45  60  75  90 

X.,(0)[d*g] 


Fig.  8.  Comparison  of  3-1)  and  2-1)  tactics. 

In  Fig.  8 the  optimal  miss  distances,  obtained  by 
a 2-1)  analysis  with  exact  (nonlinear)  kinematics 
and  the  present  3-D  study  based  on  a linearized  mo- 
del, are  compared.  The  comparison  was  carried  out 
by  a 6 degrees  of  freedom  simulation  which  used  the 
opt  i mill  control  functions  derived  by  the  respective 
studies.  The  comparison  shows  that,  for  the  set  of 
nondimens ional  parameters  chosen,  the  3-D  tactics 
have  a definite  advantage  If  the  initial  target 
azimuth  angle  Xf  > 45°. 

The  existence  of  an  optimal  maneuver  plane  en- 
ables to  use  some  results  of  the  2-1)  analysis10  and 
as  a consequence  avoids  the  solution  of  two-point 
boundary  value  problems,  which  seems  a priori  neces- 
sary if  missile  saturation  and  limited  airplane  roll- 
rate  are  considered.  By  the  way.  It  can  he  noted 
that  the  optimal  miss  distances  predicted  by  the 
three-dimensional  linearized  kinematic  model  compare 
very  well  with  results  of  complete  (nonlinear) 

6 degrees.  Such  a good  agreement  is  not  unexpected. 

It  lias  been  shown  previously10,  that  the  "bang- 
hang"  nature  of  the  optimal  mane uve r seems  to  justi- 
fy the  assumption  of  linearized  kinematics.  However 
It  is  proposed  to  distinguish  between  two  phases 
when  defending  the  choice  of  a linearized  kinematic 
model:  a)  a priori  justification  can  he  based  on 
examination  of  the  value  of  the  "dynamic  similarity 


•v  (a?)  max 

yT  = “TV, — 


(92) 


If  this  parameter  Is  sufficiently  small  use  of 
linearized  kinematics  can  be  attempted  for  the 
analysis,  b)  a posteriori , It  has  to  be  verified 
that  the  solution  of  the  linearized  kinematic  model 
does  not  predict  excessively  long  maneuvers  In  one 
direction.  If  it  does  happen,  trajectory  lineari- 
zation is  not  appropriate  for  the  specific  problem. 
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